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A pion distribution amplitude, derived from nonlocal QCD sum rules, has been 
employed to calculate J ? 7 * 7 ^ 7r (Q 2 ) using light-cone sum rules, and F n (Q 2 ) in NLO 
QCD perturbation theory. Predictions are presented for both observables and found 
to be in good agreement with the corresponding data. Calculating the hard pion 
form factor by Analytic Perturbation Theory to two-loop order, it is shown that 
the renormalization-scheme and scale-setting dependencies are diminished. 

1 Introduction 

Large-distance QCD remains an area, where the concepts of perturbation theory cannot 
be directly applied. To assess this region and make reliable predictions for hadronic 
processes, the pure perturbative treatment has to be amended by nonperturbative input. 

'Invited plenary talk presented by the first author at Hadron Structure and QCD: from Low to High 
Energies, St. Petersburg, Repino, Russia, 18-22 May 2004. 
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In a series of recent papers 012], three of us have outlined an approach, based on 
QCD sum rules with nonlocal condensates capable of providing a pion distribution 
amplitude (DA) compatible at la with the CLEO data 4 on the pion-photon transition. 
The key feature of this pion DA is that its endpoint regions x — 0, 1 (x being the parton's 
longitudinal momentum fraction) are strongly suppressed. This suppression is controlled 
by the nonlocality of the scalar quark condensate, parameterized by the average quark 
virtuality A 2 , in the vacuum, with theoretical estimates in the range (0.4 - 0.5) GeV 2 
and a preferable value of 0.4 GeV 2 extracted in [2] from the CLEO data. 

In addition, one can improve the quality of perturbatively calculable observables, 
notably the factorized hard contribution of the pion's electromagnetic form factor, by 
trading the traditional power-series perturbative expansion for a non-power-series (in an 
analytic) coupling expansion that avoids eo ipso the Landau singularity rendering all 
expressions infrared (IR) finite OE]. Suffice it here to say that this is achieved through 
the inclusion into the running coupling of a power-behaved term of nonperturbative origin 
that removes the Landau ghost leaving the ultraviolet behavior of the effective coupling 
unchanged. Crucial for making this analytic approach possible, is the generalization of 
the analytic running-coupling concept, proposed by Shirkov and Solovtsov [Hj, to the 
level of observables depending on more than one scheme scales jH|, as is, for example, 
the case for the pion form factor in fixed-order perturbation theory beyond LO or 
performing a Sudakov resummation |10| . 

Along these lines of thoughts, we describe in this contribution our recent works on 
the pion DA, summarizing the main results, and present predictions for the pion's elec- 
tromagnetic form factor carried out under the imposition of analyticity of the running 
coupling and its powers using two different procedures. It turns out that if the powers of 
the coupling have their own analytic (dispersive) images, the factorizable hard part of the 
form factor so calculated bears a minimal dependence on the scheme and scale-setting 
choice. Including also the soft contribution via local duality, this helps improving the 
quality of the prediction beyond the level of the current experimental-data accuracy. 



2 Endpoint-suppressed pion distribution amplitude 



In the context of factorization of hard exclusive processes the pion DA is a universal, 
gauge-invariant quantity defined at the twist-2 level by 
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where dxip^ (a;,/i 2 ) =].,/„■ = 131 MeV is the pion decay constant, and C(0, z) = 
Vexp[— ig^Q t a A^(y) dy^] preserves gauge invariance. ip^{x,^ 2 ) encapsulates the non- 
perturbative QCD pion structure in terms of the distribution of the longitudinal momen- 
tum fractions between its two valence partons: quark (x) and antiquark (x = 1 — x). 
Together with the DA of its first resonance, A\, it can be related to the nonlocal conden- 
sates by means of a sum rule, based on the correlator of two axial currents (see pQ). Due 
to the finiteness of the vacuum correlation length A^ 1 , the end-point regions x — > 0, 1 
are strongly suppressed and by virtue of this fact we can pQ determine quite accurately 
the first ten moments = Jq 1 (f 7T (x)(2x — 1) N dx of the pion DA and independently 

also the inverse moment (x^ 1 )^ = J Q ^(x) x~ 1 dx. Given that (£ N ) n — » rapidly 
with increasing N, the eigenfunctions decomposition 

(2) 

can be practically truncated at 04 because all higher coefficients are negligible £Q. The 
"bunch" of the pion DAs shown in Fig. 1(a), parameterized by ai and 04, turns out to 
match all moment constraints for (£, n )tt and (x -1 }^ extracted from the CLEO data. The 
optimum sample out of this "bunch" — BMS model — £Q, has at [Xq = 1 GeV a% = 0.20 and 
04 = —0.14 and is shown in Fig. 1(a). Let us close this section with a forward-looking 
statement: the BMS "bunch" pion DAs, though doubly peaked, have their endpoints 
(x — ► 0, 1) strongly suppressed — not only relative to ip^ z but even compared to <^ s , 
substantially reducing the importance of Sudakov effects. 

3 Comparison with the CLEO data on the pion-photon transition 

It was shown in ^H] at LO and later extended PJ[ to NLO of QCD perturbation theory 
[T7] that the light-cone QCD sum-rule method allows to perform all calculations in the 
l*{Q 2 )l{l 2 ) ~ ¥ 71-0 form factor for sufficiently large q 2 and analytically continue the 
results to the limit q 2 = 0, hence avoiding problems arising when a photon becomes real. 
Recently we have revised and refined this sort of data processing accounting for a 
correct ERBL ^T] evolution of the pion DA, including thresholds effects in the running 
coupling, estimating more accurately the contribution of the twist-4 contribution, and 
improving the error estimates in determining the la- and 2er-error contours. Avoiding 
here technical details, we gather the results of our analysis in Fig. 2. The predictions 
shown correspond to the following pion DA models with associated a deviations and 
designations for the form-factor predictions displayed in the right panel: ipez (■, 4cr, 



<p*(x, Mo) = &b(1 - x) 1 + a 2 (Mo) Cl' 2 {2x - 1) + a 4 (/^) C°^(2x - 1) + 



-(3/2, 




Figure 1: (a) Profiles of pion DAs normalized at \i% = 1 GeV 2 : BMS model [T]— solid 
line; CZ model [12] — dashed line; asymptotic DA — dotted line, (b) BMS "bunch" [T] in 
comparison with tp as , (pp^ (dashed line) [13], (fiuor (dash-dotted line) |14j . 



upper dashed line) BMS- "nonlocal QCD SRs bunch" (shaded rectangle), (fBMS (*, 
lcr — left panel, shaded strip — right panel) £Q; three instanton-based models, viz., ^B] 
(*, 3<7, dotted line), 2<r, dash-dotted line), and ^1 (A, 3u — only left panel); 

and the asymptotic pion DA ip aa (♦, 3er, lower dashed line). A recent transverse lattice 
result [20j (▼, 2cr) is also shown — left panel only. 

To summarize, the main results obtained in |2] are: (i) Both DAs, <p^ ^1] and (p^ z 
[T2| are disfavored by the CLEO data at 3cr and 4cr, respectively. In contrast, (fi^ MS lies 
within the la-error ellipse. Model DAs from instanton-based approaches |13M14imflll9| 
are close to but still outside the 2er region, (ii) The extracted coefficients ai and 04 are 
rather sensitive to the strong radiative corrections and the size of the twist-4 contribution, 
(hi) The value of the vacuum nonlocality extracted from the CLEO data is A 2 < 0.4 GeV 2 . 
Turning to the form-factor predictions, one observes from Fig. 2 (Right) that the BMS 
"bunch" of pion DAs is in good agreement with the CLEO data 0] but also with the 
CELLO data [2J , while the behavior of rival DAs reflects the situation shown in the left 
panel: the prediction from the CZ model overshoots the data considerably, while that 
from (p^ s — and DAs close to it — are underestimating both sets of experimental data. 

4 Electromagnetic pion form factor. Theory and phenomenology 

The crucial new elements of the calculation below are: (i) use of the BMS pion DA, (ii) 
application of two-loop Analytic Perturbation Theory (APT), and (iii) a more accurate 
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Figure 2: (a) <5 2 i ? 7 * 7 ^ 7r (Q 2 ) CLEO-data analysis in terms of error contours in the (02,04) 
plane. The line assignments are: broken line — 1<t; solid line — 2<r; dash-dotted line — 3<r. 
Various pion DAs are shown, evaluated at /x| Y = 5.76 GeV 2 after NLO ERBL evolu- 
tion.. The slanted shaded rectangle represents the nonlocal QCD sum-rule constraints on 
(02, 04) P for A 2 = 0.4 GeV 2 . (b) Light-cone sum-rule predictions for Q 2 F 1 * 1 ^(Q 2 ) 
in comparison with the CELLO (diamonds, J2I]) and the CLEO (triangles, 0]) data 
evaluated with 5 2 w _ 4 = 0.19 GeV 2 



way, based on local duality, to join the soft part with the hard form-factor contribution. 

The pion's electromagnetic form factor can be generically written as JT] Ft T (Q 2 ) = 
F Fact {Q 2 ) + F™ n - Fact (Q 2 ), where F Fact (Q 2 ) is the factorized part within pQCD and 
i^" on_Fact (Q 2 ) is the "soft" part containing subleading power-behaved (e.g., twist-4) con- 
tributions originating from nonperturbative effects. The leading-twist factorizable con- 
tribution can be expressed as a convolution in the form Fj act (Q 2 ; /x 2 ^) = f 2 ip^(x, /x|) ® 
Th(x, y, Q 2 ; ix|, /XpJ (£></?„• (j/, /xp), where /xp is the factorization scale between the long- and 
short-distance dynamics, /xr stands for the renormalization scale. The hard-scattering 
amplitude, Th(x, y, Q 2 ; /xp, /x^), describing short-distance interactions at the parton level, 
has been evaluated to NLO accuracy ([22J and references cited therein) using the termi- 
nology introduced in [H] to which we refer for details. Then, one obtains F Fact (Q 2 ; /i^) = 
F^°(Q 2 ; /x R ) + F^ LO (Q 2 ; xt|), where the LO and NLO terms read, respectively, 

^ lo (Q 2 ;a4) = «s(4)^ LO (Q 2 ) (3) 

r 1 2 

Q 2 ^°(Q 2 ) ee 8nf 2 [l + a ^ NLO (Q 2 ) + a^ NLO (0 2 )] , (4) 

F^°(Q 2 ;/ x|) = *M [^' NLO (Q 2 ;^) + ^ ND - NLO (Q 2 ;^M ax = 00)]. (5) 
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Here iVMax marks the maximal number of Gegenbauer harmonics taken into account 
and the calligraphic designation denotes quantities with their a s -dependence pulled out. 
Note that because we take into account the NLO evolution of the pion DA, the displayed 
terms contain diagonal (D) as well as (the NLO term) non-diagonal (ND) components. 
The effects of the LO DA evolution are crucial [52], while the NLO ones are relatively of 
less importance. Hence, we set here: a^' NLO — > a^' LO and a^ D ' NLO — > 0. Studying 
beyond the LO requires an optimal renormalization scheme and scale setting in order to 
minimize the influence of higher-order loop corrections (see ^9 for a fully fledged discus- 
sion). To join the hard with the soft contribution (the latter being calculated with the 
aid of local duality (LD), we have to correct the low-Q 2 behavior of the factorizable part 
to fulfill the Ward identity at Q 2 = 0, i.e., F W (Q 2 ; = F^ D {Q 2 ) + F^ act - WI (Q 2 ;/x^) 

with F^~™(Q 2 ; M |) = Q 4 / (2,s^ loop + Q 2 ' " 



if act (Q 2 ;/4) 
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Figure 3: (a) Predictions for Q 2 F 7r (Q 2 ) obtained with the BMS pion DA using standard 
pQCD within the MS scheme and adopting /x^ = Q 2 (dashed line). The solid line 
corresponds to the BLM scale setting introduced in jHI . The experimental data are taken 
from [231 (diamonds) and [23 (triangles), (b) Prediction for Q 2 F n (Q 2 ) calculated with 
the "Maximally Analytic " procedure and with the BMS "bunch" of pion DAs. 



The next step is to apply for the calculation of F^ act (Q 2 ) APT. This is done by 
employing two different analytization procedures: (i) A Maximally Analytic prescription 
0, meaning that analyticity has been imposed not only on the coupling, but also on its 
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powers, which, therefore, have their own dispersive images. This amounts to 



[j£**(Q a ;A&)] 



1 



??>m + -4 2) (^)^ NLO (Q 2 ;^) , (6) 

7T 



MaxAn 



where Sg 2 ^(/i^) is the two-loop analytic coupling and A^\/j^) the analytic version of its 
second power in two-loop order [5]. (ii) Another procedure, we call Naive Analytic, 
replaces in F^ act the strong coupling and its powers by the analytic coupling a s and its 



powers 



a s 2 (Mr) > entailing the requirement ^0] 



[^(g 2 ;4.)] NaivAn =a( 2 >(/4)^°(Q 2 ) + i [4 2) (A)] ^ nlo (Q 2 ;a4)- (7) 

The results for vs. the experimental data are displayed in Fig. 3(b) and Fig. 4. 
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Figure 4: Predictions for Q 2 F V (Q 2 ) using APT and the BMS DA in conjunction with 
the "Naive Analytic" (a) and "Maximally Analytic" (b) procedures: MS scheme and 

/Up = Q 2 (dashed line); BLM (dotted line); BLM (solid line); ay-scheme (dash-dotted 
line). The single solid line in panel (b) shows the prediction for the soft form- factor part; 
below this, the corresponding hard contributions are displayed. 



5 Summary and conclusions 

The BMS pion DAs pQ successfully pass the comparison with the CLEO data 0] at the la 
level, as highlighted in Fig. 2 (conforming also with the CELLO data Employing 2- 

loop APT — naive and maximal — we have calculated the hard part of the electromagnetic 
pion form factor within various renormalization schemes and using different scale settings. 
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Joining the hard part with the soft one on the basis of local duality, we have derived 
predictions that reproduce the available data rather well, especially using the "Maximally 
Analytic" procedure (Fig. 3(b)). Moreover, we found that this procedure minimizes the 
influence of scheme and scale-setting ambiguities on the form- factor predictions. 
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